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I. INTRODUCTION 

In recent years the study of non-nuclear vulnerability and 
survivability of aircraft flying in a hostile arena has become increas- 
ingly important. Increased aircraft sophistification has brought about 
increased complexity of aircraft systems. Better surface-to-air 
defenses, high weapon system costs and decreasing defense budgets have 
made survivability analysis and design absolutely necessary. 

Aircraft fuel cells are very susceptible to damage because they 
have the largest surface area and volume of all vulnerable aircraft 
components. Ballistic threats have caused extensive fuel cell structural 
damage which can result in fuel starvation, fires or explosion. Major 
fuel cell structural damage cannot be tolerated in the case of "wet 
wing" aircraft for obvious reasons. It is for these reasons that the 
study of the failure mechanisms of fuel cells must be understood so 
that survivable designs may be evolved. 

Projectiles which penetrate fluid filled cells cause damage many 
times more severe than that incurred by impact with an empty cell. This 
interaction of fluid filled cells with projectiles is called hydraulic 
ram. The exact nature and extent of damage can be related to projectile 
energy loss in each of five phases of hydraulic ram. These phases are 
usually identified by the dominant energy transfer mechanism presented. 

The first phase involves projectile penetration through the entry 
wall material. During this entry phase projectile energy is lost by 
"punching" a hole in the entry wall material. If the projectile is not 
tumbling and strikes normal to the entry wall the damage size is 
approximately equal to the maximum cross sectional area of the threat. 
Tumbling and incident strikes can produce larger damage sizes and 
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correspondingly higher penetration energy loss. 

If sufficient projectile energy is available the threat protrudes 
from the back side of the entry wall and subsequent fluid displacement 
produces a strong hemispherical shock wave which propagates into the 
tank fluid. This shock phase of hydraulic ram generates intense 
pressure loadings on the entry wall but because of spatial attentuation 
these loadings are localized to the immediate area surrounding the 
impact point. If the wall is thin., and the threat kinetic energy is 
large, or if during the penetration phase the entry wall has been 
sever ly cracked, shock phase pressure loadings can cause extensive 
damage to the entry wall area. 

As the projectile travels through the tank fluid a cylindrical 
cavity is formed. The cavity is a result of a radial velocity imparted 
to the fluid by the projectile motion. Inside this cavity are trapped 
air, which entered during penetration, and fluid vapor. Projectile 
form drag dissipates projectile energy in this drag phase of hydraulic 
ram. 

If the fuel cell has a depth which is less than the distance re- 
quired for the drag phase to dissipate all projectile energy (for 
example the case of a projectile penetrating a "wet wing" from the lower 
surface) an exit phase of hydraulic ram exists. The exit phase in- 
volves exit wall pressure loadings due to projectile wall proximity and 
due to fluid motion toward the wall that has been established by the 
earlier phases. These loadings usually cause extensive, large scale 
damage to the exit wall. 

The last phase is produced when the drag phase cavity growth is 
halted by fluid pressures due to the tank boundaries . The cylindrical 
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cavity changes to a spherical shape and collapses. Collapse compresses 

the trapped cavity vapors to a high pressure resulting in a re-expansion 

of the cavity. Expansion and collapse is repeated several times before 

the stored energy is dissipated. The resulting fluid pressures can cause 

large side wall loadings which can further damage the cell. These 

oscillatory cavity phase pressures can destroy large depth cells. The 

cavity phase is also responsible for pumping tank fluid out damaged 

« 

areas of small depth tanks. This release of fuel, atomized by wall 
cracks and petaling, creates an extreme fire hazard. 

All phases of hydraulic ram must be considered when designing "ram 
proof" fuel cells. Fuel cell damage and the dominant hydraulic ram 
failure mechanism are dependent upon the aircraft mission (which deter- 
mines the path, energy level and size of the most probable threat) and 
the geometry of the fuel cell (size, location, and construction). For 
example, if the tank wall construction is "thin" and the threat energy 
is "high" the shock phase of hydraulic ram can cause catastrophic damage 
of the entry wall. Tanks with "small" depth are very susceptible to 
damage produced during the exit phase of hydraulic ram while the cavity 
phase contributes only minor additional structural damage. Fuselage 
tanks may have dimensions sufficient to eliminate the exit wall phase. 

The primary ram damage component is the cavity phase in this case. 
Efficient fuel cell protection must be tailored to the dominant ram 
components. Since ram damage is configuration dependent extensive 
vulnerability testing will be required to substantiate the defeat of a 
given threat. 

This report presents the results of experimental studies conducted 
at the Naval Postgraduate School to understand the basic phenomena 
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associated with the first three phases of hydraulic 



ram. 
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II. PENETRATION PHASE OF HYDRAULIC RAM 



The assessment of ram damage to any fuel cell can be related to 
the energy dissipated in each of the five phases of hydraulic ram. For 
this reason it is necessary to be able to predict the projectile energy 
loss during entry wall penetration. In general this energy loss is a 
function of wall parameters and projectile initial energy, obliquity, 
material properties, tumbling and shape. For the purposes of supporting 
NPS experimental studies of the shock and drag phases of hydraulic ram, 
it was necessary to be able to predict the penetration energy require- 
ments of untumbled, copper- jacketed, small arms lead projectiles with 
various masses and nose shapes when impacting normal to 7075 -T 6 aluminum 
tank entry walls of various thicknesses. 

Dunn's^ theory for hypervelocity ballistic impact problems assumes 
that "there exists in each unconstrained material body at constant 
uniform temperature a dynamic yield strength, a , that is a constant of 

y 

the material." This dynamic yield strength when multiplied by the 
resisting cross sectional area of the projectile gives an estimate of 
the instantaneous retarding force during penetration. The projectile 
equation of motion 

x + -£ A (x) = 0 (II-1) 

m 

may then be solved to yield the time history of penetration depth and 
projectile velocity assuming that the dynamic yield stress of the pro- 
jectile is greater than that of the entry wall. The work done during 
deceleration yields an estimate of the kinetic energy loss during 
penetration. For the retarding force postulated, the work done is 
equal to the product of the dynamic yield stress and the volume of 
entry wall material swept out by the penetration. The theory, therefore, 
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predicts that the penetration energy loss is independent of projectile 
shape. The minimum initial projectile velocity which will penetrate 
the entry wall is called the ballistic limit and is given by: 



V = 
o 



2 a At 

y 

m 

\ 



(II-2) 



where A is the damage area and t is the entry wall thickness. 

The exit velocity V g of the projectile penetrating a given wall 
thickness is given by: 




The corresponding percentage energy loss is given by: 
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(II-4) 



Experimental studies conducted by Dunn suggest that the correct value 



of a is approximately fifteen times the shear modulus of the entry wall. 

y 



These theoretical predictions were compared to a series of experi- 
2 

ments in which projectile shape, impact velocity, and mass were varied. 
The tank entry wall was simulated by a 7075-T6 Aluminum plate of thick- 
ness 0.05, 0.09, or 0.l6 inches. Three projectile nose shapes ranging 
from hemispherical to a 1.6 fineness ratio ogive were tested. Figure 
II-l shows the components of the ballistic range used for these 
experiments. The rifle mount shown in Figure II -2 was adjustable in 
azimuth and elevation and was designed to hold rifles up to and including 
30 caliber. The impact velocity was measured by using three chronographs 
screens mounted in front of the test plate. Each screen had a five volt 
D.C. source across it. When the screen was broken by the passage of the 
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bullet the circuitry shown in Figure II-3 shaped the pulse and started 
or stopped a Monsanto 101B lMH z counter. Time interval between screens 
was converted to projectile average velocity to determine the initial 
velocity. A similar pair of chronograph screens were mounted behind 
the test plate to measure the exit velocity. Figure II-4 shows a typi- 
cal chronograph station. 

Early in the test phase it was experimentally determined that the 
size and method of supporting the test plate had no significant effect 
on the exit velocity. It is expected therefore, that the exit velocity 
for hypervelocity impact would be insensitive to wall-fluid foundation 
effects that would be present in any fluid filled fuel cell. Several 
rounds were fired through a single 18 -inch square plate and the 
resulting exit velocities were identical even if the point of impact 
was as close as two times the damage diameter away from any previous 
damage. All further shots used a three by four inch test plate to 
simulate the fuel tank entry wall. Figure II- 5 shows the experimental 
apparatus . 

Spark shadowgraphs of the penetration were also taken. The delay 
circuitry shown in Figure II-3 was used to trigger the spark source 
just after the penetration. Figure II-6 shows the optical setup. The 
light from the spark source was directed to a film plate by a concave 
mirror on the opposite wall. The collimated light then passed behind 
the test plate and onto the film plate. These shadowgraphs showed 
that, except for minor distortion of the nose tip, that these copper- 
jacketed lead projectiles show no appreciable deformation during the 
penetration. Figure II-7 shows a typical shadowgraph. The yield 
strength of copper is approximately 0.39 X 10^ psi as compared to 
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0.25 X 10 psi for Aluminum. The theory assumes that the yield 
strength of the projectile is significantly larger than that of the 
entry wall. Since these data show that projectile shape was un- 
changed it can be concluded that the assumption only need be made that 
the yield strengths of the projectile is equal to or greater than that 
of the entry wall for hypervelocity impact problems. 

Figure II-8 shows a comparison of experimental impact and exit 
velocities with the theory for a 45 grain 22.2 caliber projectile 
impacting plates of various thicknesses. Projectile nose shape was 
found to have no effect on the results ; therefore , no differentiation 
between shapes has been made on this figure. It is obvious that 
equations II-2 and II-3 predict very accurately the exit velocity for 
hyperballistic impact. Figure II-9 shows the resulting penetration 
energy comparison. Good agreement with the prediction equation II-4 
is obtained. Additional shots with bullet masses as large as 63 
grains have shown similar excellent agreement with theory. 

The independence of penetration energy loss with projectile shape 
makes it possible to extend these results to many cases of interest to 
the fuel cell designer. It could be expected, for example, that the 
penetration energy of warhead fragments impacting normal to any inter- 
vening structure and the fuel cell wall can be predicted accurately 
using this simple theory. 
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FIGURE II-l. Ballistic Range Components. 
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FIGJTPB II-3* Adjustable Time Delay Unit Circuitry 
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FIGURE 11-^. Ballistic Range Chronograph Station 
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FIGURE IIt 5. Test Plate and Chronograph Station 
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FIGURE II -6. Spark Source and Mirror 





FIGURE II-7- 
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IMPACT VELOCITY (ft/sec) 

FIGURE II-9* Impact Velocity vs. Percent Kinetic Energy Loss for Aluminum 
Plates of Three Thicknesses 



III. Shock Phase of Hydraulic Ram 

In order to predict fuel cell damage due to the shock phase of 

hydraulic ram it is necessary to relate the change in projectile energy 

3 

to wall pressure loading. The results of several investigators have 
shown that protrusion of the projectile into the tank fluid produces a 
hemispherical shock wave which is centered at the impact point. Pro- 
jectile momentum and path after entry have no observable effect on the 
shape of this shock wave as it propagates through the tank fluid. In 
fact, even impacts with obliquity show a hemispherical shock front 
centered at the impact point. It may be concluded from these results 
that the shock characteristics may be predicted by assuming that the 
shock is the result of a point energy source release. 

To predict fuel cell internal pressures created by such an event 

the computer code in the Appendix was developed based on the theory of 

5 

Yurkovitch . The theory assumes that the tank walls are rigid and that 
the shock radius is proportional to the 0.8 power of time. The one- 
dimensional Rankine-Hugoniot equations are coupled to the fluid state 
equation to determine shock front conditions. The equation of state is 
based upon the Tait equation for isothermal compressible liquids when 
fluid changes are produced by an adiabatic and isentropic process. The 
assumed state equation was: 



P = 



P C 
o o 

n 



l_ K o 






(III-l) 



where p Q is the density and C q the speed of sound in the undisturbed 
fluid and n is a constant for a given fluid (n = 7*0 for water, n = 10.6 
for aviation fuel. 

Nondimensionalization of the variables was accomplished by the 
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following scheme: 



5 = 57*7 



f(5,M ) = p ^ ?t) 

p o 

0(S> M s ) - 4t*} (III - 2) 

s 

p o K s <*> 



s(5.«.) - 



#*> 

where R is the shock radius, u is the particle velocity, p is the 
s 

pressure and e is the specific internal energy of the fluid. By using 
the strong shock assumptions that the shock front pressure and energy 
are much greater than the undisturbed conditions the equations of 
motion reduce to: 



1 - \ (1-^) 
f l = ^l 

&L = |f L (1 - | ) (III-3) 

f = -i-*- (Y n -l) 

nM 

s 
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Solving for the conditions at the shock front 




s 



0 - (Y n -l) 

nM 

s 



(III- 1 *) 



f l-*L 

and assuming values of Y, the set III-4 can be solved for M , f and 

-L S 1 

0^ respectively. The shock front fluid parameters are now known as a 
function of the shock Mach number, M , but their distribution behind the 
front was to be determined. 

The distribution of the flow parameters behind the shock is based 
on the assumption that the density profile is given by 

q( M_) 

Y = Y 1 ? (III-5) 

The exponent q(M ) is determined by equating the total mass enclosed 
s 

by the shock 

3 1 2 

P 0 R s 2 tt Io (III “ 6) 

to the original mass occupied by the enclosed volume 



2riR^p 
s r o 

3 

with the result 

Substitution of (ill- 5) into (III-8) yields 



(III-7) 



(III -8) 
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(HI-9) 



q(M g ) =» 3 (\ -1) 

the nondimens ional particle velocity 0 was calculated by substituting 
the density profile into the continuity equation 



w n + * n + 2 ^ m (in-io) 

^ ^ s 

2 

where p is the decay coefficient given by R R /R . The resulting 

s s s 

particle velocity was 

0 = (III-11) 

The known particle velocity and density distributions were then 
substituted into the momentum equation 



f - f i - jrfe (s'^hr - e Wi - M s tr> 



iA - - 

to solve for the pressure distribution. The decay coefficient was 
determined by assuming that the shock radius is given by 



(III-12) 



R = Ct 
s 



0.8 



(III-13) 



The shock Mach number was determined from energy considerations . 
The energy imparted to the fluid in the shock phase, E q , is given by 



R 



E o " $o (e - e o + P 4 ™" dr 



(III-14) 



Using the normal strong shock assumptions and the previously derived 
density and particle velocity distribution equation III-lU becomes 

E o • 01 + 01 (XXX-15) 

nM 
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By choosing M the dependent varibles 0 and Y.. are known and equation 

S XX 

o 

III-15 can be integrated to yield 2E /R° versus Mach number. Given the 

o' s 

shock phase energy release E q it is then possible to determine the shock 
Mach number as a function of shock radius. There is, however, an anomaly 
at a shock radius of zero. To avoid a finite energy in zero volume 
reference 6 suggests that the initial radius be given as 



R 



o 




2 np C 2 
r o o 



1/3 



(111-16) 
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A series of experiments were conducted to substantiate the 
theoretical propagation speed of the shock wave in a water filled tank. 
Figure III-l shows the ballistic range components used for these tests. 
Projectile impact kinetic energy was measured using three chronograph 
screens which triggered counters. The time interval between screens 
was converted to average velocity and then to impact energy. A fourth 
chronograph screen was mounted ahead of the test tank entry wall. This 
screen provided a reference pulse for timing events during the shock 
phase. The circuity in Figure II-3 was used to trigger a pair of spark 
sources at preset times after the reference pulse. In this manner, a 
pair of spark shadowgraphs of the events during the shock phase was 
obtained. Figure III-2 shows the optical arrangement of the shadowgraph 
apparatus . 

Figure III-3 shows the test tank and stand. The fuel cell was 
approximately cubical with 17 inch inside dimensions. The front and 
back wall support plates were machined one inch thick aluminum. These 
plates had a 15 inch diameter circular cut-out. The simulated fuel 
cell entry and exit walls were mounted inside these plates with 16 
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equally spaced bolts. The entry wall was a .090 inch thick 7075-T6 
aluminum sheet with a prepunched entry hole as shown in Figure III-4. 

The top of the tank was open and tape was used to prevent water leakage 
before impact. The side and bottom walls of the tank were one inch 
thick ELexiglass on which a reference gird was scribed. 

The experiments were conducted at three initial projectile 
energy levels. Figures III-5 and III-6 show the 22 caliber projectile 
shapes and energy parameters. By varying the preset time delay shadow- 
graphs of the shock radius versus time were taken. Figure III-7 shows 
typical shadowgraphs . The times shown in this figure are after wall 
impact. The hemispherical shock centered at the impact point is readily 
apparent. The projectile is obscured by light refraction caused by the 

large density gradients produced in the stagnation region of the flow. 

3 

McMillian has shown, using X-ray shadowgrams, that the illuminated cusp 
can be used to determine the location of the projectile. Also apparent 
in the later time shadowgraph is the cylindrical cavity caused by flow 
separation from the projectile surface. It should also be noted that 
these shadowgraphs show that the shock strength (as indicated by the 
band width) is not uniform. The strength is a maximum directly ahead 
of the projectile with a weakening as the wall is approached. This 
effect was probably due to projectile forward motion and entry wall 
motion in responce to the shock pressure loadings. 

Figure III-8 through III-10 show a comparison of the experimental 
results with the theory. The data for an impact energy of 1600 in-lb 
shows excellent agreement with the theory assumtion that the shock 
radius varies with the 0.8 power of time. The higher energy level data 
show, however, that the shock radius varies as 0.9 power of time. This 
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may be due to the differences in projectile shape. If the projectile is 

considered to be a piston in the early stages of the shock formation, 

the rate of change of fluid volume displaced by the piston determines 

the momentum exchange between the piston and the fluid, hence the 

shock strength. The low energy projectile had a nose shape that was 

more blunt than the higher energy projectiles. It is reasonable to 

conclude therefore, that the greater rate of change of fluid volume is 

displaced by the blunt projectile producing a stronger shock. The 0.8 

4 

power law has been experimentally verified by Stepka for the case of 
spherical projectiles. The present data shows that this law may be ex- 
tended to the case of blunt projectiles or possibly warhead fragments. 
The data also indicates that a 0.9 power law more closely predicts the 
shock propagation characteristics for projectiles with nose shapes that 
approximate an ogive. 

Also shown in Figure III-10 (X symbols) are data for the case of 
projectile penetration through a solid 0.043 inch thick 6061-T6 aluminum 
entry wall. These limited results show that the shock propagation speed 
is relatively insensitive to the entry wall penetration process. The 
spark shadowgraphs in Figure III-ll show, however, that the solid wall 
shock phase wave system is considerably more complex than the prepunched 
case. Entry wall motion in response to the penetration phase produces 
a weak wave system with the shape of a truncated cone which intersects 
the main shock tangentially. The base dimension of this weaker shock 
was found to be approximately equal to the product of the speed of 
sound in aluminum and time. The primary shook radius, however, was 
found to be approximately equal to that of the prepunched wall case. 

The resulting wave interactions completely obscure the projectile in 
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this case and presumably alter the wall pressure loadings during the 
shock phase. Figures III-12 through l4 show the damage resulting from 
this type of test. These figures show that although insufficient energy 
was available to cause entry wall catastropic failure, considerable 
plastic deformation and buckling occur ed. Since the exit wall exhibited 
no damage, it is reasonable to assume that the entry wall damage was 
due to shock phase loadings with minor contributions due to loadings 
produced by the latter phases of hydraulic ram. 

Implicit in the Yurkovitch solution for shock phase pressure 
loadings is a knowledge of the energy input E q . An analytical expression 
for this energy as a function of projectile, wall, and fluid characteris- 
tics is not available. In addition, the motion of the thin fuel cell 
walls in response to these impressed pressures will reduce the fluid 
pressure in the vicinity of the wall. An estimate of these effects is 
one of the ultimate goals of the present NPS hydraulic ram research 
effort. 

If, in the design effort to build shock phase survivable fuel cells, 
the most conservative approach is assumed (all of the projectile energy 
is lost and the walls are rigid) the solution of Yurkovitch can yield 
entry wall pressure loadings until such a time that the shock wave re- 
flects from any other tank wall. Figures III-15 through III-21 show 
such calculated pressure profiles for the three energy levels of the 
previous experimental data. Figures III-15, l6, 18 and 20 show the pre- 
dicted pressure signals that would exist at a constant radial position 
verses time. These figures show that the shock strength is quickly 
attentuated spatially. It is expected, therefore, that shock phase 
damage would be limited to the immediate area surrounding the impact 
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point. Figure III-22 shows the shock radius and time at which the pre- 
dicted shock wave Mach number is almost sonic (M = 1.05). The theory 
assumes no counter pressure (p Q « p). For this reason it would be 
expected that sonic shock speed would actually be reached at an earlier 
time and smaller radius than shown. However, this figure provides a 
conservative estimate of the possible extent of damaging shock phase 
pressure loadings and of the time scale involved in the shock phase 
phenomenum. These theoretical predictions for the time at which the 
shock is acoustic were 22 micro-seconds at 1.71 inches for the 1600 
in-lb energy level, 36 micro-seconds at 2.31 inches for the 7^93 in-lb 
energy and 44 micro-seconds at 5.4 inches for the 12323 in-lb energy 
level. 
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figure III- 1 . Ballistic Range Components 





FIGURE III-2. Schematic of Shadowgraph Apparatus 

(Down-range View) 
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FIGURE III-3. 



Shadowbox, 



Tank and 



Stand Installation 
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FIGURE III-4. Front Wall Test Specimen, 
Fre-punched Hole 
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igure IIJ-5. Projectile Shapes (L.to R . .22 
caliber, m=5.75 x I0 -3 lb, .222 caliber, m= 6.43 
x|0“3|b, .222 caliber, m = 7.85 x I0" 3 lb) 



Caliber 


.22 Long Rifle 


.22 2 Remington 


.222 Remington 


Impact Velocity, Vo 
(ft /sec) 


1200 


2500 


2900 


Projectile Mass, m 
<lb m ) 


5,75xl0“ 3 


6.43xl0 -3 


7.85 x I0” 3 


Impact Energy, E 0 
(in-lb) 


1600 


7493 


12323 


Experimental Drag 
Coefficient, Cq 


1.20 


0.45 


i 

0.50 


Frontal Area of 
Projectile After 
Traversing Fluid, A 
(in 2 ) 


0.122 


0.098 


0.165 


Equivalent Flat 
Plate Area, f (=C n A) 

M 


0.146 


0.044 


0.083 



igure IIT-6. Projectile Parameters 
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Figure III-7- Shadowgraphs After Projectile Impact Through 

a Pre-punc'ned Front Wall 
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Theoretical 
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III-8. shock Radius vs. Time, = 1600 in-lb 



Theoretical 
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Figure ITT-ST. Shock Radius vs. Time, E_. = 74 93 in-lb 



Solid Front Wall 
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Figure m-10. shock Radius vs. Time, E„ = 12323 in-lb 




Figure III-ll. Shock Waves Resulting from 0.22 Caliber 
Projectile Impact with a Solid Entry 
Wall 
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Figure III-12 . Solid Front Wall £ f ter Test (6061-T6 Aluminum, 
.043 in. thick, E Q = 12323 in-lb) 




Figure III-13. Typical Damage to Solid Front Wall at Point of 
Impact (Front view, 6061-T6 Aluminum, .043 
thick, E c = 12323 in-lb) 
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Figure III-14., Damaoe to folid Front V'all (Mounted on 







Pressure - P xio 5 (psi) 




Figure III-15. Pressure vs. Radius from Impact - 
Complete Profile, E Q = 1600 in-lb 
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Pressure - P xlO “ 4 (psi) 




Figure UI-16. Pressure vs. Radius from Impact, 

= 1600 in-lb 
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figure III-17. Pressure vs. Time for Three Radii, E n = 1600 in-lb 
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’igure IH-l8. Pressure vs. Radius from Impact, 
E Q = 7493 in-lb 
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Figure III-19. Pressure vs. Time for Three Radii, E = 74 93 in-lb 
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Figure HI-20. Pressure vs. Radius from Impact, 

E = 12323 in-lb 
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Time - t (micro seconds) 

Figure III -22. Radius and Time for Acoustic Shock 
Speed in Water Verses Energy Input 
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IV. DRAG PHASE OF HYDRAULIC RAM 



After protrusion from the entry wall the projectile transfers its 
remaining energy to the fluid by the production of the shock phase pre- 
viously discussed and by viscous dissipation as it traverses the tank 
fluid. Figure IV-1 presents a qualitative view of this projectile 
energy decay. The shock phase does not, in most cases of interest, 
dissipate all of the projectile energy. The projectile subsequently 
travels through the tank until the remaining energy is dissipated or 
until the exit wall is reached. During this latter drag phase of 
hydraulic ram a cylindrical cavity is formed as a result of radial 
velocity imparted to the fluid by the projectile motion. Inside this 
cavity are trapped air which entered during the entry wall penetration 
and fluid vapors. 

During the drag phase projectile spin and angle of attack combine 
to produce Magnus sideforces which perturb the projectile from a straight 
path. Furthermore, large stagnation pressures in the flow cause projec- 
tile deformation and in some cases the projectile does not remain intact. 
It is also likely that the energy imparted to the fluid by viscous 

dissipation during the shock phase would modify the pressure distribu- 

Q 

tion behind the shock. Fuhs has suggested that a characteristic time 
separating the drag and shock phases of hydraulic ram be identified as 
the point in time where projectile stagnation pressure just exceeds the 
shock front pressure. This characteristic time was approximately 30 
micro-seconds for the 1600 in-lb energy level projectile tested. Inter- 
actions between the shock phase and drag phase pressure loadings would 
be largest for times close to this characteristic time. 

These factors make an analytical estimation of drag phase pressure 
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loading exceedingly difficult. Stepka has suggested that during the 
drag phase the projectile energy decay may he approximated by solutions 
of a simple equation of motion 

D - - m g - C D \ p^A (IV-1) 

In general projectile deformation would require both drag coefficient 
and area variation with time. By assuming constant projectile area, 
drag coefficient, and mass equation IV-1 may be integrated to yield the 
projectile velocity decay as a function of distance traveled through the 
fluid 




-Bx 



e 



(IV-2) 



where B = pC.pA/2m and is the initial velocity of the drag phase. 
Projectile position as a function of time is given by 



x » \ in (BV Q t + 1) (IV-3) 

From this crude analysis an estimate of the rate of energy transfer to 
the fluid may be made from: 



E_ ^ 1 

E 0 [1 + BV 0 t] 2 



(IV-4) 



7 

A series of experiments were conducted to determine the drag phase 
energy dissipation. Using the experimental setup discussed in section 
III spark shadowgraphs of the drag phase were taken for three initial 
impact energy levels. Figure IV-2 shows typical results. The cylin- 
drical cavity formation is clearly shown in these shadowgraphs. The tip 
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of the illuminated cusp marks the position of the projectile nose. From 
measurements of projectile position versus time the energy decay of the 
drag phase may he determined. Figure IV-3 through TV-5 show a comparison 
of the experimental and theoretical positions for the three energy 
levels tested. It was hoped that this data would yield an estimate of 
the energy loss during the shock phase, hence the unknown initial drag 
phase velocity V Q defined in Figure IV-1. Data scatter, however, was 
too great to make such an estimate and was subsequently assumed to 
be equal to the impact velocity. Data scatter was attributed to varia- 
tions, from shot to shot, in projectile path and deformation. It was 
concluded that trajectory parameters hence drag phase energy decay is 
very sensitive to projectile attitude and deformation during fluid 
traverse. 

Approximately 51% of the 22 caliber projectile impact velocity was 
lost within 1.5 inches of the front wall and within 135 micro-seconds. 

The 22.2 caliber projectile with an initial energy of 7^93 in-lb lost 
only 17% of its impact velocity at 1.5 inches while those with an energy 
of 12,323 in-lb lost 25%. The differences in these results are attri- 
buted to differences in projectile deformation rates. The lead 22 
caliber blunt projectile deforms more rapidly that the copper jacketed 
22.2 caliber rounds. It would be expected that this round would lose a 
larger percentage of the initial energy at a constant penetration depth 
because of the resulting larger equivalent flat plate area (f = C^A). 

The percentage energy lost at a constant radius was also a function of 
the initial energy level for projectiles of the same construction. In- 
creased initial energy produces larger viscous dissipation rates and 
larger deformation rates because of higher stagnation pressure loadings. 
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This effect is shown clearly by examining typical projectile shapes 
after fluid traverse as shown in Figure IV-6. The projectiles studied 
are known to by dynamically stable when the longitudinal axis is perpen- 
dicular to the flight path and hence they do not tumble. Stagnation 
pressures generated by the drag phase flow are large enough to cause 
large plastic deformations of these lead and copper jacketed projectiles. 
Since drag coefficient and area are changing during the drag phase it is 
more convenient to use an equivalent flat plate area to match experimen- 
tal data to the predictions of equation IV-4. These results indicate 
that the exponential velocity decay can be used for a reasonable estimate 
of the drag phase energy loss if an equivalent flat plate area of the 
round in question can be determined experimentally. These data also 
indicate that the drag phase energy dissipation rate is very sensitive 
to the material and aerodynamic properties of the threat. 
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.Figure 17-0.. Qualitative View of Projectile Energy and Velocity Decay in the Fluid 





Figure IV-2. Cavity Growth Resulting from a 
Prepunched Hole Impact 
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Theoretical 
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Figure TV-3. Projectile Distance Traveled in Fluid vs. Time, E Q = 1600 in-lb 
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Figure IV-4. Projectile Distance in Fluid vs. Time, E 0 = 7493 in-lb 
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Figure IV-5. Projectile Distance in Fluid vs. Time, E = 12323 in-lb 




Figure IV-6. Typical Deformation of Projectiles After 
Coming to Rest in the Fluid Filled Tank - 
Pre-punched Hole Results (L. to R. : .22 

caliber, Ei = 1600 in-lb; .222 caliber, 

E = 7493 in- lb; .222 caliber, E 0 = 12323 
in-lb) 
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V. Recommendations for Further Study 

Hydraulic ram studies will continue in FY75 and the following list 
outlines the future work: 

1. Information similar to that presented in this report will 
be collected using the higher energy levels of 30 _ caliber 
projectiles. 

2. The test cell will be instrumented with pressure trans- 
ducers to determine experimental pressure profiles during 
the shock phase of hydraulic ram for cells with rigid walls. 

A comparison of these results with predictions will allow 
an assessment of the Yurkovitch model and could provide a 
method for determining the shock phase energy loss E q . 

3. Fuel cell pressure loadings will be determined for cells 
with thin entry walls. A comparison with the results in 
item 2 should result in a determination of the fluid-wall 
motion interaction problem. 

4. The test cell entry walls will be instrumented with 
strain gages to determine wall stresses during hydraulic ram. 

These data will be compared to theoretical predictions based 
upon both the experimental and the theoretical pressure 
loadings. The theoretical predictions will be obtained from 
the digital computer code SATANS 10 and BR-1HR. 
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VIII. APPENDIX 



COMPUTER PROGRAM 



THIS PROGRAM CALCULATES THE FOLLOWING 



PART 

1. 

2. 

PART 

1 



I 



RADIUS VERSES TIME 
MACH NUMBER VERSES TIME 



******************************************** ** *********** 
* 

4 
* 

* 

* 

* 

* 

* 

4 
* 

* 

* 

* 

* 

* 

* 

* 

* 

********************************************************* 



SHOCK 
SHOCK 
II 

PRESSURE VERSE TIME AND RADIUS 
FROM IMPACT BEHIND THE ShCCK 
FCR A GIVEN FLUID AND IMPACT ENERGY, EO. 
ASSUMES: SHOCK RADIUS IS PROPORTIONAL TO 
C.8 POWER UNTIL ACOUSTIC, A POWER LAW DEN 
BEHIND THE SHOCK, A STRONG SHOCK, ADIA6AT 
THE SHOCK, THE UNDISTURBED FLUID PRESSURE 
TIALLY ZERO, SPECIFIC INTERNAL ENERGY IN 
FLUID IS ESSENTIALLY ZERO. 



THIS PROGRAM 
TIME TO THE 
SITY PROFILE 
IC ACROSS 
IS ESSEN- 
LNCI STURBEC 



SYMBOL DESIGNATION 

EO= IMPACT ENERGY OF PROJECTILE ( IN-LE ) 

CO = SPEED OF SOUND FOR FLUIC (FPS) 

RHOO= DENSITY IN UNDISTURBED FLUID ( LE/FT3 ) 

DC= INCREMENT OF NON-DIMENSIONAL RADIUS, DELTA 
ZHETA 

DELT= INCREMENT OF TIME (MICRC-SEC) 

ET= ENERGY/VOL. = 2*E0/(R**3) 

FI= NON-DIMENSIONAL PRESSURE AT THE SHOCK VS. 
SHCCK MACH NO. 

DODM= DERIVATIVE OF NON-DIMENSIONAL PARTICLE 
VELOCITY W.R.T. SHOCK MACH NO. VS ShCCK 
MACH NO. 

VS= SHCCK MACH NC. 

X= PERCENT EO LOST AT ENTRY WALL 



DIMENSION ET(20) , PSI 1 ( 20 ) , FI ( 20 ) ,VS(20) , CCDM(20) 

1 , RXMS (20),CRS(20),ARMS(20)»DT(20)»T(20) ,RE(20), 

1R V ( 20) , ROOD (20) , XM S (20 ) , RS ( 20) , R PS I (20 ) , R F (20 ) 
DIMENSION TT (300) ,XM(300) , R( 300 ) ,F2( 300) , PSI2 (300) , 
1DCCM2( 300 ) , Q (300 ) ,C(100) ,F(2 ,20 ) ,P(2 ,20) 

DIMENSION C0NRS(30C) ,SHGCV(300) ,TCF0R(3 00) 

CIMENSICN FOR (20 ) 

CC 500 J=1 , 20 
FOR ( J ) = 0.0 
5CC CONTINUE 
N = 1 8 

READ (5,100) EO»CO,RHOO,X,DC,DELT 
READ (5 ,101) (ET( I ) .1 = 1, N ) 

READ (5,101)(PSI1(I) , 1 = 1 , N ) 

RE AC ( 5 , 101 ) ( FI ( I ) , 1 = 1, N ) 

READ (5,101) ( VS( I ) , 1=1, N ) 

REAC ( 5 , 101 ) ( DOOM ( I ) , I = 1 , N) 

WRITE (6,200 ) EO , CO , RhCO ,X,DC,DELT 
DC 15 1=1, N 

15 WRITE (6,201) ET ( I ) , PS 1 1 ( I ) , F 1 ( I ) , VS( I ) ,DCCM ( I ) 

PART I 
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PERCENTAGE uF INITIAL ENERGY DEPOSITED ANC INITIAL 
SHCCK RADIUS 

EG=X*EO 

R C=(E0*32-. 2*144. / ( RHOO*CO*CO*2 .*3.1416) ) **0.33333 
INITIAL VALUE OF ENERGY/VOL. 

ETO=2.*EC/(RO**3.0) 

INVERT INPUT ARRAYS FOR USE WITH INTERPOLATION 
ROUTINE, PIF2 

DC 12 1*1, N 
J=N-I+1 
RE ( I ) = ET ( J ) 

R PS I ( I )=PSI1(J) 

R F ( I ) = F 1( J ) 

RV(I)*VS( J) 

12 RCCD ( I ) =DGDM ( J ) 

INTERPOLATE FOR INITIAL VALUES AT SHOCK 

VS0=PIF2(ET0,RE,N,RV) 

F10=PIF2(VS0,RV,N, RF) 

PSI10 = P IF2(VS0,RV,N,RPSI) 

OCDMiO = P I F2 ( VSO , RV , N , RDCC ) 

COMPUTE INITIAL TIME FOR DISTANCE RO 

TC=R0/(CC*VS 0*12.0) 

ASSIGN SUBSEQUENT VALUES OF SHCCK MACH NC . AT WHICH 
SHOCK RADIUS WILL BE COMPUTED 

DC 10 1*1, N 

IF ( VSO-VS ( I ) ) 1,2,2 

1 GC TO 10 

2 M=N-I+1 
XMS(1)=VS0 
GC TO 11 

1C CONTINUE 

COMPUTE SUBSEQUENT VALUES OF SHCCK RADIUS 

11 RS(1)=R0 

DO 20 J=2,M 

XMS ( J ) = VS.( I) 

RSI J) = (2.*E0 /(ET (I)) ) **0 . 3333 

1 = 1 + 1 

20 CCNT INUE 
K = M-1 

COMPUTE DELTA TIME CORRESPONDING TO SHOCK MACH NO. 
ANC SHOCK RADIUS 

CC 30 1*1, K 

DR S 1 1 ) = RS(I + 1)-RS(I) 

ARMS! I ) =2-/1 XMS ( I ) +XMS I 1 + 1 ) ) 

3 0 DT(I)=DPS(I)*ARMS( I)/(C0*12.0) 

T 1 1 ) = T0 
DC 40 J=2, M 

40 T ( J ) = T ( J— 1 ) +DT I J-l ) 

WRITE (6,102) 

WRITE (6,103) EO 
WRITE (6,110) 

PRINT SHGCK MACH NC. AND SHCCK RACIUS VS TIME 
DC 50 1*1, M 

5 C hRITE (6.111) T ( I ) , XMS ( I ) , RS ( I ) 

100 FORMAT ( 6E10.0 ) 
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I C 1 FORMAT (8E10.0) 

102 FCRMAT (1H,//,39X, 'ENERGY ' , 15X , 'FLUID* ) 

1C3 FCRMAT (29X,1E20.7) 

11C FORMAT (1H,//,39X 'TIME' »17X,'MACH' »17X , ' RACIUS' ) 
111 FCRMAT ( 29X » 3E20 . 7 1 
2CC FCRMAT (//, 6E20.7) 

2C1 FORMAT (15X,5E20.7) 



PART II 



INITIALIZE VARIABLES AND ASSIGN VALUES USED WITH 
NCN-DI M ENS IONAL PRESSURE DISTRIBUTION FUNCTION 

8EB=T(M)*1.E06+1.05 
L = EBB 
TT ( 1) =0.0 
XM ( 1)=XMS (1) 

R ( 1 ) = R0 
F 2 ( 1 ) =F 10 
PS 12 ( 1 ) =PS 1 10 
DCDM2(1 )=DCDM0 
Q(1)=3.*<PSI10-1.0) 

E=-0 . 23 

1=1. 0/DC +1.05 

cm=o.o 

WRITE (6,301) 

WRITE (6,302) TT ( 1 ) , XM ( 1 ) , R ( 1 ) 

WRITE (6,303) 

DC 60 J=2,N 

6C C(J)=C( J-1)+DC 

COMPUTE NON-DIMENSICNAL AND CIMENSICNAL PRESSURE FOR 
INITIAL ShCCK RADIUS. 

DC 70 J=1.I 

F( 1,J)=F2( 1)+PSI2( 1)*(C( J)**(0( l)+2.0)-l.C)*(F2(l)/ 
1PSI2(1)-B*(F2(1)+XM(1)*DGDM2(1 ) ) ) / (Q( 1 ) +2 .0 ) 

P( 1,J)=F(1,J)*RHOO*CO*CO*XM(1)*XM(1)/ (32 .2*144. ) 

PUNCH 91,R(1),C(J),P(1,J) 

7C CONTINUE 
K= 1 
JJ = 1 
SLM=0. 0 

CALL FORCE(FOR,C,P,R,I,JJ,K) 

DC 501 J=1,I 
SLM=SUM+FQR ( J ) 

WRITE (6,304) C(J),P(1,J),F(1,J),FCR(J) 

5C1 CCNTINUE 

TCFOR ( 1 )=SUM 

ITERATION TO SOLVE PRESSURE, SHCCK RADIUS, SHCCK MACH 
NO. VS TIME 

DC 80 K=2,L 

TT(K)=TT(K-1)+DELT 

IF (TT(K)-TO) 21,21,22 

21 XM ( K) =XM (J. ) 

R(K)=R( 1) 

F2(K)=F2(1) 

PSI2(K)=PSI2(1) 

DC DM2 ( K ) =D0DM2 ( 1 ) 

0 ( K ) = Q ( 1 ) 

CC 55 J=1,I 

55 P(2,J)=P(1,J) 

PUNCH 92 , R ( K) 

GC TO 80 

22 X=TT ( K ) 

INTERPOLATION OF INPUT DATA AND PART I DATA FCR 
VARIABLES IN PRESSURE FUNCTION 
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XM K) = P I F2 ( X , 1,M,XMS) 

R(K)=PIF2(X,T,M,RS) 

Y = XM ( K ) 

F2(K)=PIF2(Y,RV,N,PF) 

FSI2(K)=PI F2(Y,RV,N,RPSI) 

DCDM2(K)=PIF2 (Y . P V , N , RDCD) 

Q(K)=3.*(PSi2 (K)-l.C) 

RPITE (6,301) 

RRITE (6,302) TT(K),XM(K),R(K) 

RR ITE (6,303) 

DC SO J=1,I 

F(2,J)=F2(K)+PSI2(K)*(C(J)**(Q(K)+2.0)-1.0)*(F2(K)/ 
1PSI2(K)*-B*(F2(K)+XM(K ) * DO DM 2 ( K ) ))/(Q(K)+2.0) 
P(2,J)=F(2,J)*RH0C*C0*C0*XM(K)*XM(K)/ (32.2*144. ) 

PUNCH 91, R C K ) ,C( J) ,P( 2, J) 

CONTINUE 
J J-2 
S0M=0.0 

CAIL FCRCE(FGP,C,P,R, I , J J , K ) 

DC 502 J = 1,I 
SLN=SUM+FGP( J ) 



5o2 

8C 



91 

92 
2 C 1 

302 

303 
3C4 
3C5 



RPITfdg^CH) C(J) ,F(2,J),F(2,J) ,FCR(J) 



C 

C 

C 



9 

3C 

10 

11 

12 

13 

14 

15 

16 

17 

18 



TCFCR ( K ) =SLM 
CONTINUE 
RPITE ( 6 ,305 ) 

CALL CQNVER(R,XM,TCFGR,TT,L) 

FGRMAT (3E15.5) 

FCPMAT (1E15.5) 

FGRMAT (1H»// ,36X, 'TINE' ,17X,'MACF' ,17X, 'PACIUS' ) 
FCPMAT ( 29X, 3E20 . 7 ) 

FCPMAT ( 1H , // ,39X, ' P /RS ' , 17X , • P • , 17X , • F • , 1SX, 'FORCE* ) 
FCPMAT(29X»4E20.7) 

FCPMAT ( lh,//,36X, • TIME' ,17X, 'VELC ' ,17X, 'RACLIS', 17X, *T 
1GFGRCE ' ) 

END 

FUNCTION PIF2 ( X , X L 1ST , N , FL 1ST ) 

SECOND GRDER INTERPOLATION 

DIMENSION XL I ST (100), FLIST (100) 

ELIF (P,Q,R,S,T) = ( (Q-P)*(S-T)/(R-Q)+S) 

IF (X-XLIST(N)) 2,1,1 
I = N- 4 1 
GC TO 5 

I F ( X-XL 1ST ( 1 ) ) 4,4,6 
I = 1 
K = 1 
GC TO 3C 
K = 2 

DO 8 I = 1 , N 

IF (X— XL 1ST < I ) ) 9,9,8 

CONTINUE 

I =N 

I = 1-1 

ELI FI = BLIF(X,XLIST(I),XLIST(I + 1),FLIST(I),FLIST(I + 1) 
1) 

IF (K-l) 11,11,12 
PIF2 = BLIF1 
RETURN 

I F ( ( I +2 )-N ) 13,13,16 
IF ((I-l)-l) 15,14,14 

I F ( ABS ( XL 1ST ( 1-1 )-X)-ABS( XL 1ST (1 + 2 )-X)) 16,15,15 
L = 1+2 
GG TO 17 
L = 1-1 

BLIF2 = BLIF (X, XL 1ST ( I ) , XL 1ST ( L ) , FLIST ( I ) , FL 1ST ( L ) ) 
PIF2 = oLIF (X,XLIST( 1+1) ,XLIST(L) ,BLIF1,BLIF2) 

RETURN 

ENC 
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SLBROUT INE FORCE ( FOR,C ,P,R, I , JJ ,K) 
******************************************* ************** 



* 

* 

* 

* 

* 

a 

<r 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



THIS SUBROUTINE COMPUTES THE FORCE ON THE ENTRY HALL 
DUE TO SHOCK PHASE FOR A FRONT HALL WITH A 
PRE-PUNCH.ED HOLE , I . E. PROJECTILE COES NOT IMPINGE CN 
FRONT WALL UPON ENTRY. COMPUTES FORCE TC EE ZERO FCR* 
PRESSURE FIELD OVER THE PRE-PUNCHEC HOLE CP WHERE THE* 
PRESSURE GOES NEGATIVE. THE PROGRAM ACCOMPLISHES THIS* 
BY A STAIGHT LINE APPROXIMATION EETWEEN PRESSURE * 

DATA POINTS. * 

RHRS= NON-DIMENSIONAL RADIUS OF PRE-PUNCHEC HOLE * 

S= SLOPE OF STRAIGHT LINE BETWEEN PRESSURE CATA * 

POINTS * 

RX = THE INTERCEPT OF R/ RS , NON-C I M ENS I ON A L RACIUS, * 
WHERE PRESSURE IS ZERO * 

FOR (N ) = FORCE (LB ) FOR AN ANNULUS WHOSE OUTSIDE 
IS C ( N) AND INSIDE RADIUS IS C(N-l) 



RADIUS* 
* 

* 

************************** **************** *************** 



CIMENSICN FOR (20) , C ( 1 00 ) t P( 2 1 20 ) ,R(300) 

RFRS = .3 75/ R ( K ) 

I F ( RHRS-1.. ) 100 » 25 » 25 
IOC IF(JJ-l) 3,1,2 

1 CC 10 N=2, I 

S=(P( 1 ,N)-P( 1 , N— 1 ) )*4. 

IF(P(1,N-D) 3,4,4 

3 IF ( P ( 1 ,N) ) 5,5,6 

5 FOR ( N ) = C.O 
C-C TC 10 

6 RX=(S*C(N-1)-P(1,N-1) )/S 
IF(RHRS .GT.C(N) ) GO TO 43 
IF(RX.GT.RHRS) GO TC 60 
Ph=S*(RHRS-C(N-l) )+P( 1,N-1) 

FCR(N)=( .5*(P(1,N)-PF)+PH)*3.1416*(C(N)*C(N)-RHRS 
1*RHRS)*R(K)*R(K) 

C-C TO 1C 

60 FCR(N)=P(1,N)*.5*3.1416*(C<N)*C(N)-RX*RX)*R(K)*R(K) 
GO TO 10 

4 IF(RHRS-C(N-1) ) 40,40,41 

AC FCR(N)=(.5*(P(1»N)-P( 1 , N-l ) ) +P ( 1 , N-l ) ) *3 . 1 A 16* ( C (N ) 
1*C(N)-C(N-1)*C(N-1) )*R(K)*R(K) 

GO TO 10 

A 1 IF (RHRS-C ( N) ) 42,43,43 
42 PF=S*(RHRS-C(N-1) )+P( 1,N— 1) 

FCR(N)=(.5*(P(1»N)-Fh)+PH)*3.1416*(C(N)*C(N)-RHRS 

1*RHRS)*R(K)#R(K) 

GC TO 10 
A3 FOR(N) =0.0 
1C CONTINUE 
RETURN 

2 DC 20 N *2 , I 
S=(P(2,N)-P(2,N-1) )*4. 

I F ( P ( 2 , N— .1 ) ) 13,14,14 

12 IF(P(2,N)) 15,15,16 

15 FCR (N) =0.0 
GC TO 20 

16 RX=(S*C(N-1)— P(2,N-1) )/S 
IF(RHRS.GT.C(N) ) GC TO 143 
IF ( RX. GT.RHRS ) GO TO 160 

PH = S*(RhRS-C(N-l ) )+P(2,N-l) 

FCR(N)=(.5*(P(2,N)-PF)+PH)*3.1416*(C(N)*C(N J-PHRS 
1*RHRS)*R(K)*R(K) 

GC TO 20 

16C FCR(N) = P(2,N)*.5*3.l4l6*(C(N)*C(N )-RX*RX ) *R (K)*R(K) 
GC TO 20 

1A I F ( RHRS-C ( N-l ) ) 140,140,141 
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140 FCR(N)=(.5*(P(2,NJ~P(2,N-1) )+P(2,N-l) )*3.1416*(C(N) 
1*C(N)-C(N-1)*C(N-1 ))*R(K)*R(K) 

C-C TO 2C 

1 1 IF (RHRS-C(N) ) 142,143,143 

142 PH=S*(RHRS-C(N-1) )+P(2,N-l) 
FCR(N)=(.5*(P(2,N)-PH)+PH)*3.i41fc*(C(N)*C(N)-RHRS 

1*RHRS)*P(K)*R(K) 

GC TO 2 G 

143 FCR (N) =G.O 
2C CONTINUE 
25 RETURN 

ENC 

SUBROUTINE CONVER ( R , XM , TOFOR , TT , L ) 

THIS SUBROUTINE CONVERTS VELOCITY AND RACIUS 
TC KM/SEC AND CM RESPECTIVELY 

CIMENSICN R(3G0) , X M 300 ), TOFOR ( 300 ), CONRS ( 3 CO ), SHOC V ( 
13CC ) , TT (300 
DC 510 11=1, L 
CCNRS(II)=R(II)*2.54 
S HOC V ( I I ) = XM ( II ) *4. 9/3. 281 

WRITE (6,520) TT( II ) » SHOCV (II ) , CONRS (I I ) , TCFCRdI ) 

52C FORMAT (29X,4E20.7) 

510 CONTINUE 
RETURN 
END 
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